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Abstract. Using results derived in Berestycki et al. [T] we correct the 
celebrated formulae of Hagan et al. [5] H]. We derive explicitly the cor- 
rect zero order term in the expansion of the implied volatility in time to 
maturity. The new term is consistent as /3 ^ 1. Furthermore, numerical 
simulations show that it reduces or eliminates known pathologies of the 
earlier formula. 



We discuss hei^ the celebrated formulae of Hagan and Woodward [5] and 
Hagan et al. [4J, point out some inconsistencies and compare them with the 
results of Berestycki et al. [1]. This leads to a new corrected version of the 
original formula derived in [3]. This is a short technical paper and the reader 
is assumed to be familiar with the SABR model and asymptotic expansions of 
the implied volatility surface. We refer to Gatheral [_3j Chp.7] for background. 

Consider the following model for the underlying (we take r = which should 
be interpreted as St representing the forward value process) 

dSt = atS'^dWl, So = s, 

dat = uatdWt, do = a, where d{W^,W'^)t = pdt, 

(3 G (0, 1), IpI < 1, which is known as the SABR model. The model ([1]) is 
written in the pricing measure under which are correlated Brownian 

motions. We take v > and note that putting u = we obtain a local 
volatility model. We refer interchangeably to strike X or to x = ln{s/K). The 
implied volatility is denoted I{x,t), i.e. a European call with strike K and 
time to maturity r has the same price under the SABR model ^ and under 
the Black-Scholes model with volatility I{x,t). 

Consider now the Taylor expansion of the implied volatility surface I{x,t) 
in time to maturity r 

lix, t) = I\x) f 1 + I\x)t) + 0(r2). (2) 
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Table 1. Comparison of term in [4J and [Ij. 



Hagan et al. [3] were first to obtain explicit expressions for 1° and They 
used perturbation theory (combined with impressive intuition) . More recently, 
Berestycki et al. jLlj treated the subject in a rigorous analytical manner. They 
proved in particular that I^{x) = lim^^o t) is well defined and character- 
ized it via a solution to eikonal equation (or via Varadhan's signed geodesic 
distance). Direct, but somewhat tedious, computations allow us to compare 
their expression for with that of [1] and of [5j (case of = 0). The findings 
are summarised in Table 1. 

The formula for I^ix) in the local volatility case v = Q agrees in both works 
and appears to be correct. We note however that it differs from the formula in 
Hagan and Woodward [3]. Numerically the difference is small but persistent. 
Unfortunately, it is the formula in [S] which is often used, as in Bourgade and 
Croissant P]. 

The formula for I^{x) when /? = 1 is the same in |4j and in py. However the 
two papers diffeiS when (3 <1 and we believe it is the formula of [1] which is 
correclo and should be used. This is natural as authors in pQ present an ana- 
lytic treatment of the subject and characterize theoretically. Furthermore, 
the form of their formula we derived here coincides with the general formula 
obtained by Labordere [6l Eq. (5.1)] (taken for the case of 0-SABR). We give 
now two arguments in favor of the formula of [Ij: theoretical and numerical. 

From the theoretical point of view, the formula of Hagan et al. [1] has an 
internal flaw: it is inconsistent as (3 ^ 1. Indeed, I^{x) in the formula of [3j 
does not converge to the known value for I^[x)p=i as /? — > 10 In contrast, as 

V s^~^ — K^^^ vx 

z = z{(5) = > — , when ^ ^ 1, 

a 1 — p a 

we have convergence using the formula of [T] . 



^Note there is a typo in (6.5) in [T] which we corrected here. 

"^Strictly speaking Hagan et al. |4] do not claim to obtain the Taylor expansion in maturity, 
however their formula (A. 65) is of the form ([2]). 

^However, when one uses the simplified formula (2.17a), which is (A. 69c), instead of the 
general formula (A. 65), the convergence issue (magically) disappears. 
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The formula of Hagan et al. [4J is known to produce wrong prices in region 
of small strikes for large maturities. It assigns a negative price to a structure 
with a positive payoff, or equivalently it implies a negative probability density 
for the stock price process in some region. If one uses the formula we derive 
here the problem either appears in yet lower strikes or disappears completely. 
To see this, we price the following structure in interest rates: long 3% put and 
short 1% put at three times the notional, whose payoff is a triangle, zero at 
zero and at 3% and with a peak at 1%. We denote it T(l). More generally, 
T{K) is the following: long {K + 2)% put and short iC% put at {K + 2)/K 
times notional. The payoff oiT{K) is a triangle, zero at zero and at (if + 2)%, 
with a peak at K%. 

Let and I%{x) be the formulae for from ^ and pLj respectively. 

The implied volatility Ih{x) in [U Eq. (A. 65)] is given via 

lH{x) = I%{l + lUx)r), 

1 {P-lf a" 1 pvaP 2-3p^ ^ (3) 

'^^^^ 24 4(,^)(i-W2+ 24 

and we define the implied volatility Ib{x) = + /]^(x)r), as we did 

not derive an explicit formulae for l]^ from [1]. An explicit formula for if , 
albeit complicated, is derived in Labordere [B], who also argues that I}j IS a 
valid approximation. Figures 1 — 2 compare the prices of the structure 




T{K) derived using Ih and Ib- The anomalies do not occur, or occur at lower 
strikes, if we price using the implied volatility Is instead of Ih- 
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Figure 2. Implied volatilities and T{K) prices for maturity r = 20 
p = -37%, /? = 0.6, ATM = 4.25%, u = 24.5%, s = 8.01% 
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